Fermion localization on degenerate and critical branes 



R. A. C. Correa* A. de Souza DutraJ M. B. Hott* 



UNESP Univ Estadual Paulista - Campus de Guaratingueta - DFQ. 
Av. Dr. Ariberto Pereira Cunha, 333 
12516-410 Guaratingueta SP Brasil 

March 28, 2011 



Abstract 

In this work we analyze the localization of fermions on degenerate and critical Bloch 
branes. This is done directly on physical coordinates, in constrast to some works that 
has been using conformal coordinates. We find the range of coupling constants of the 
interaction of fermions with the scalar fields that allow us to have normalizable fermion 
zero-mode localized on the brane on both, critical and degenerate Bloch branes. In the 
case of critical branes our results agree with those found in [Class. Quantum Grav. 
27 (2010) 185001]. The results on fermion localization on degenerate Bloch branes 
are new. We also propose a coupling of fermions to the scalar fields which leads to 
localization of massless fermion on both sides of a double-brane. 



1 Introduction 

The idea that the Universe we live in can be realized by a static 3-domain wall (3-brane) 
immersed in a (4,l)-dimensional world has opened a pathway for the localization of matter [lj 
and gauge bosons [2] in worlds with large extra dimensions without resorting to mechanisms 
of compactification [3]. Large extra dimensions has also provided mechanisms to solve the 
hierarchy of interactions problem [4]-[7j as well as the cosmological constant problem [8]. 
In [9] it was shown that the effective gravitational potential between two particles recovers 
the Newtonian behavior, since one has localization of gravitons on a thin brane in five- 
dimensional space-time with an warped geometry and the cosmological constant is related 
to the brane tension. Later the localization of matter (spin-zero, spin- 1/2 and spin-3/2) in 
the Randall- Sundrun framework was shown to be possible, under certain conditions over the 
brane tension [TU]. It is important to emphasize that the localization of fermions on thin 
branes is, in fact, provided by an ad hoc soliton and the mechanism for localization follows 
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straightforwardly a la Jackiw and Rebbi [TT]. As a matter of fact, it is such a soliton that 
provides the domain-wall and the fermion localization in the scenario proposed by Rubakov 
and Shaposhnikov pQ. 

By introducing a non-linear model with a set of scalar fields, in 5-dimensional space- 
time with warped geometry, one has a set of coupled non-linear differential equations whose 
minimum energy solutions are thick branes and self-consistent warp factors. The thick 
branes also separate the space in two patches characterized by a peculiar warp factor whose 
asymptotic behavior is an anti-de Sitter (AdS5) space [12]- [19], as in the Randall-Sundrun 
framework (thin brane), but without singularities. The stability of such solutions is very 
difficult to be proven due to the intricate differential equations one has to deal with [T7] . 
Notwithstanding, one still has localization of gravitons on thick branes, as has been shown 
in [18]. In fact, this program has been developed as a generalization of the domain wall 
universe, by taking into account the stabilization of gravity fluctuations via domain-walls in 
supergravity theories [2TJ . 

Localization of matter on thick branes has been illustrated by using several different 
non-linear models for scalar fields coupled with gravity [22]- [21]. One of those models is 
the Bloch branes model [22], which comprises two interacting scalar fields whose classical 
solutions are the Bogomolnyi-Prasad-Sommerfield (BPS) and the warp factor can also be 
obtained as solution of a first-order differential equation. Moreover, the model exhibits a 
richer structure due to the variety of kinks (solitons) it comprises [2S]-[2H], leading to what 
has been called degenerate and critical Bloch branes |29j, with a self-consistent warp factor 
and localization of gravitons. A natural track that has been followed is the localization of 
matter in such a variety of branes. 

In [30] the localization of massless fermions has been studied together with an analysis 
of resonant fermion modes with a specific coupling of fermions to the classical BPS config- 
urations of Bloch branes and in [31] it is studied the localization of massless fermions on 
critical Bloch branes. It has been explicitly pointed out [21] that critical Bloch branes arise 
as a critical limit of degenerate Bloch branes due to a running, but limited, constant of the 
integration of the orbit equation relating the classical configurations for the scalar fields. 
Thus the localization of massless fermions on degenerate and critical branes is an important 
question to be analyzed. As a matter of fact, a Bloch brane may be seen as a thick brane 
that evolutes to a thicker one, namely, a degenerate Bloch brane which, in its turn, splits 
into two branes whose separation becomes large as the degeneracy parameter approaches 
the critical value, at which the critical brane is triggered. Such a picture resembles the de- 
scription of first order phase transitions which was also used in the context of brane worlds 
to describe brane splitting [32]. In fact, as a counterpoint, in the model we use here, the 
free energy itself does not depend on the temperature; instead, this dependence is implicit 
in the degeneracy parameter. 

The objective of the present work is to analyze the fate of massless fermions trapped on 
a split brane. In section II we give a brief review of the model to be used, together with the 
consistent branes and warp factors solutions. In section III we deal with the localization of 
massless fermions on the critical and degenerate branes and the IV section is devoted to the 
conclusions and further remarks. 
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2 A brief review of a model and its degenerate and 
critical Bloch branes 

The action in five-dimensional gravity coupled to two interacting real scalar fields, can be 
represented by 



S = J d 4 xdr y/\g\ 
where g = det (g a b) and 



~R + \{d^(j> + dpx&x) ~ V(<P, X) 



(1) 



ds 2 = g ab dx a dx b = e 2A{r \^dx^dx u - dr 2 , (a, b = 0, 4), (2) 

where r is the extra dimension, t]^ u the usual Minkowski metric in the four space-time 
dimensions and e 2A(r ) is the so-called warp factor. 

If the potential V(4>, x) can De written in terms of a superpotential as 



/ dW(<f>,x) \ 2 | f dW^xV 2 ' 
V dcj) ) \ dx 



^(0, X ) 2 , 



(3) 



By substituting the superpotential, and using the orbit equation, we obtain 



A(x) 

Mx) 



(A ^ 2/x), 



'2\a 2 \ 1 

«i + ( ^— I Mx) 



(3/i + Aci) 2 1 2 A , , 1 



( ln(x) -2 



(4) 



(A = 2/i), (5) 



9/i y w 18/i 

where «o and ai are arbitrary integration constants, which are chosen to be A(r = 0) = 0. 

It has also been found in (2HI EH] that the classical solutions for Co < —2a and A = fi are 
given by 



2a 2 



<t>. 



(i) 



{\/ c l ~ 4 fl2 ) cosh(2uar) — . 
a ^a/cq — 4a 2 \ sinh(2/iar) 



'Cq — 4a 2 J cosh(2«ar) — c 
The corresponding warp factor is expressed as 



(6) 



(7) 



3 2A(r) 



N 



2a 2 



^a/cq — 4a 2 ) cosh(2/iar) — c 



4a 2 /9 



X 



x exp 



2a 2 


Cq ± 4a 2 — c ^1 


/c 2 - 4a 2 ) 


cosh(2a/ir) 
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- 4« 2 ) 


cosh(2a/ir) — Co 


2 



(8) 



where N is chosen such that e 2A ^ = 1, for plotting convenience. 

On the other hand, for A = 4/x and Cq < l/(16a 2 ) the solutions can be written as 



2a 



!> (2) ( 

DBW^ 



— 16c a 2 ) cosh(4/iar) + 1 

a (a/1 — 16c a 2 ) sinh(4«ar) 
(Vl — 16coa 2 ) cosh(4/iar) + 1 



with the warp factor 



N 



2a 



— 16c a 2 ) cosh(4/iar) + 1 



16a 2 /9 



X 



x exp 



4a 2 

~9~ 



1 + 8a 2 c + — 16a 2 c ) cosh[4/xar] 
(1 + (VI - 16a 2 c ) cosh[4/xar]) 2 



(9) 
(10) 



(11) 



The set of solutions above, was baptized by Dutra and Hott [29] as degenerate Bloch 
walls (DBW). 

Furthermore, an interesting class of analytical solutions, named as critical Bloch walls 
(CBW), was shown to exist when the constant of integration equals a critical value. For 
A = \i and Cq = —2a, one has the set of solutions for the scalar fields 



XcBw( r ) — 
0Gw( r ) = 

which leads to the following warp factor 



[1 ± tanh(/xar)] 



-— [tanh[/iar) =)= 1] 



(12) 
(13) 



D 2A(r) 



■ a >, 2a 2 /9 

N\ - [1 ±tanh(/xar)]| x 



exp( — [1 — tanh 2 (/iar]). 
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For A = Afi and Cq = l/(16a 2 ) the solutions for the fields are given by 



(14) 



XcBw( r ) 
4> { CBw( r ) 



V2 a 



cosh (fiar) ± sinh (fiar) 
A/cosh(2uar) 



— [±1 — tanh(2a/ir)] 



(15) 
(16) 



and the warp factor is 



e 



2A(r) = N 



2 fl 2 e ±2^iar 1 8a2 / 9 



cosh(2/iar) 
2a 2 e ±2 ^ ar 



X 



x cxp 



1 + 



9 cosh(2/iar) [ 4cosh(2yuar) 



(17) 



In Fig.l are shown profiles of the warp factor in the case of DBW with A = 4/i for some 
values of the constant of integration cq and in Fig. 2 the behavior of the warp fact in the 
case CBW and A = 4fi is shown for some values of the parameter a. 

We would like to warn the reader that the expressions for the warp factor and the Figures 

2 and 3 presented in [29] do not agree with those presented here. The correct expressions 
and Figures are those presented here. Despite this mistake, the conclusions of the work [2S] 
are not wrong, except for the fact that we had attributed the appearance of two peaks in 
the warp factor, presented in figure 2 of that work, as a sign of the formation of two domain 
walls. As a matter of fact, two- kink solutions can be seen from the expressions ([7]) and 
({TO]) for values of cq close to the critical value. In Fig 3. we illustrate the two-kink solution 
with A = 4/x for two values of the constant of integration Cq. In brane cosmology scenario 
we are interested in, the formation of two branes leads to an almost flat space-time region 
between them, that is, the warp factor does vary appreciably between the two branes, as 
can be depicted from Figure 1. The critical brane could be seen as two branes, but infinitely 
separated from each other, such that one has a complete wetting. 

The two-kink solutions of this model were recently investigated to discuss the phe- 
nomenon of brane splitting by means of an effective model with only one scalar field [33] . In 
fact, that effective model was built based on the very same model explored here . The con- 
stant of integration Co plays the role of a coupling constant in the effective potential obtained 
in [33], such that one could think of the effective potential as a free energy that describes 
a first-order phase transition characterized by the emergence of a growing wet (disordered) 
phase in between two ordered phases [32]. The branes are the domain walls separating the 
disordered domain from the ordered ones, as can be seen from the behavior of the energy 
density of the matter fields. 

3 Localization of fermions 

In this section, we study the localization of massive fermions on the degenerate and critical 
Bloch branes [29]. For this, we consider a Dirac spinor field coupled with the scalar fields by 
a general Yukawa coupling. Thus, the action we are going to work is given by 

S1/2 = / d 5 x^=g {mr a D a y - r^F(0, , (is) 

consequently the equation of motion is 

[zT a D a ~r ] F(<f ) ,x)]^ = 0, (19) 



5 



where F(<f>, x) is a functional of the classical configurations which are solutions of the equa- 
tions (??). The equation of motion for the fermion can be rewritten as 

[zr^ + ^r 4 J D 4 -^F(0, x )]^ = o. (20) 

The relations between the warped-space gamma matrices ({r a ,r 6 } = 2g ab ), with g ab defined 
in (j2J), and the Minkowskian ones ({7^,7^} = 2r/" y ) can be realized as follows 

r M = e-^V and r 4 = -^7 5 - ( 21 ) 

Moreover, we have the following expression for the covariant derivative 

D a = (d a + u a ) = d a + KjI ~ b r s r- b , (22) 

where a and b, denote the local Lorentz indices. Thus the spin connection b is given by 



ul 1 = l -E b *{d a E$ - d b E b ) - l -E b \d a Et - d b El) - l -E c *& \d c E p - q - d p E c 9 )E*. (23) 
In the above definition E% is the vielbein, and the non-vanishing components of u a are 

= ± e A (d r A) lfll5 . (24) 
Then the equation of motion for the fermion is given by 

{iVfy + e A{r) l5 [d r + 2d r A{r)\ - rje^F^, r) = 0. (25) 

From now on we make use the general chiral decomposition 

\j?(x, r) = } u ipLn(x)a Ln (r) + ^ 'ipRn(x)a Rn (r), (26) 

n n 

with 7 5 V , l„(x) = — ipLn(x) an d l br >pRn{.x) = ipRn(x). Furthermore, we assume that ipLn(x) 
and ipR n (x) satisfy the 4-dimensional massive Dirac equations 



iYd^Ln{x) = m n ifj Rn (x), 



(27) 



i^d^Rnix) = m n ip Ln (x). 



(2f 



Thus, applying the chiral decomposition ( 12"B"]) in equation f[2"5"j) and using equations ( |2"7|) 
and PS|) . we arrive at the equations below 



d dA(r) , ' 

3; +2 -sr -**"(*•*) 

d _oL4(r 



dr dr 



+ r)F((f>,x) 



-m n e A{r) a Ln (r), 



otLn{r) = m n e A{x) a Rn (r) 



(29) 
(30) 
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for the r-dependent parts of the spinor ty(x,r). In order to find reliable results concerning 
the fermion localization in the brane we make use of the following ortonormalization relations 
for the r-dependent parts of the spinor: 

30 P OO 

e 3A(r) a Ln {r)a Lm {r)dr = / e 3A{r) a Rn {r)a Rm {r)dr = 5 nm , (31) 

oo J — oo 

oo 

e 3A{r) a Ln (r)a Rm (r)dr = 0. (32) 

■oo 

By redefining the r-dependent parts of the spinor as 

a Rn (r) = e- 2A ^L Rn (r), a Ln (r) = e~ 2A ^L Ln {r), (33) 



we are able to get rid of the second-term in the left-hand sides of the equations ( 1291) and 
(1301). Thus we obtain 



L Rn (r) = -m n e~ A ^L Ln (r), (34) 

L Ln (r) = m n e- A ^L Rn (r). (35) 

We notice that the above equations are equivalent to the equations for the components 
of a spinor describing a massless fermion in 1+1 dimensions subject to a mixing of scalar 
and vector potentials. The time-independent equation for a fermion under such potentials 
can be written as Hip{r) = Eip(r), with the Dirac Hamiltonian given (in natural units) by 
H = a 2 p + V s (r)ai + V v (r), where p = —id/dr is the momentum operator, o\ and <7 2 are the 
two non-diagonal Pauli matrices, V s (r) = —r)F((f), x) (note that F((p, x) is a function of r) is 
the scalar potential and V v {r) = m n e~ A ^ is the vector potential. In this analogy, one can say 
that L^nir) and L Rn (r) play the role of the upper and lower components for the fermion zero- 
mode in 1+1 dimensions. We have mentioned that the fermion in 1+1 dimensions is massless, 
but in fact, the scalar potential can be thought as a position-dependent mass. As one knows, 
many examples of such systems were already solved in the literature [31], particularly when 
the scalar potential is proportional to the vector potential, namely V s (r) = SV v (r), that 
allow for fermion bound states. The vector potential can be attractive for fermions whilst 
it is repulsive for antifermions and vice-versa, that is, one can have pair-production. On 
the other hand, if S > 1 and the vector potential is attractive for anti-fermions, the mixing 
of such potentials supports bound states. This could be explained due to an increasing of 
the threshold for the pair-production provided by V s (r), since it contributes as a variable 
mass for the fermion and the energy provided by the electric field does not reach two times 
the effective mass of the fermion. In the brane world scenario we are considering here, 
massive (m n ^ 0) as well as massless (m n = 0) fermions might be localized inside the brane, 
depending on the shape and strength of the coupling with the scalar field r]F((p, x), although, 
the issue of localization of massive modes in the brane is a very difficult as discussed below. 

Concerning the localization of massless fermions in the brane, one finds that the r- 
dependent coefficients appearing in the chiral decomposition fl30|) are written as 



— -r]F((f),x) 



— +r)F{(f>,x) 
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/r 
F{r')dr'\ 

/r 
F{r')dr'\ , (36) 

where we have defined F{r) = F(<f)(r),x{r)). In general, one resorts to an analytical func- 
tion F (4>(r) , x( r )) i n order to have localized massless fermions in the brane and those 
localized states have invariably a well defined chirality, which depends on the behavior of 
jr-s> oo p^ r i^ r anc j on sign of 77, since the normalization condition for the r-dependent 
parts of the massless fermion are given by 

poo poc 

\N Ro \ 2 e~ A ^ +2 ^ r F ^ dr ' dr = \N L f e -^)-j,fW dr = i. (37) 



An interesting approach of this, can be found in the work by Slatyer and Volkas [24J. In 
general, if a R (r) is normalizable, (r) is not and vice-versa, then the non-normalizable 
contribution equals to zero. This could explain why we observe neutrinos with only one 
chirality in our universe. The fact that only one of the chiralities is normalizable can be 
seen from the effective Schrodinger equation that can be derived for Lrq{t) and Llq{t) 
which are zero-modes of effective potentials that are supersymmetric partners of each other, 
and the normalization of both components would imply into the supersymmetry breaking. 
One can also check that the factor e~ A ^ in the integrands above is also important for the 
normalization of the spinor, although it is not decisive for the resolution of the chirality, 
since it amounts to the same weight in the normalization for both chiralities. 

For m n ^ we observe that the left and right components can be decoupled. It can be 
shown that L Ln {r) obeys the following second-order differential equation 

L'L + A'L' Ln + [rjA'F + r,F' - V 2 F 2 + m 2 n e- 2A }L Ln = 0, (38) 
where the prime denotes derivative with respect to r. With the redefinition 

L Ln (r)=e- A ^ 2 f Ln (r). (39) 
one finds that fLn{ r ) obeys a time- independent Schrodinger equation 

-f'L + U^ f h(r) = 0, (40) 
with eigenvalue equals zero and the effective potential 

U? ff (r) = 7] 2 F 2 - V F' - rjA'F + (1/4) A' 2 + (1/2) A" - m 2 e" 2A(r) . (41) 

For the right- component one finds 

- f'L + Uf ff f Rn (r) = 0, (42) 

with 

U? ff {r) = rfF 2 + r]F' + V A'F + (1/4)A' 2 + (1/2) A" - m 2 n e~ 2A{r) . (43) 
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We have reduced the problem of localization of massive fermions on a brane into a Sturm- 
Liouville problem, which agrees with the results found in [22J and [21] , without resorting to 
the transformation of variable as done in others papers. 

One can note that there is a symmetry relating the effective potentials, namely U^f — 
^f/| an d that the eigenstates fimir) an d fnn(r) are subject to distinct effective poten- 
tials. In order to have localization of a specific massive mode m n , each one of the chiralities 
must be a bound state of its corresponding potential with energy equal to zero. Due to the 
above symmetry and the dependence on the mass in the expressions (141 p and (|4"5|) one has 
a specific potential for each massive mode and specific chirality. In general, the equations 
(flOj) and ( H2|) are very difficult to be analytically solved, but there are two particular and 
unphysical cases, namely rjF = m n e~ A ^ r ' and rjF = —m n e~ A<yr \ for which one can find exact 
expressions for aR n (r) and a^ n (r), namely an n (r) = «L n (r) ~ exp[— 2A(r)}, which are not 
normalizable. 

Now, in order to compare our results with a previous one [31], we address the issue of 
possible localization of massless fermions on the brane by setting the general coupling 

F((j), X) = ui<p + u 2 x + w 3 <fix- (44) 

to the fermions, where the u^s are constant parameters to be determined such that (r) 
or (r) is normalizable. We have analyzed the scenarios of both degenerate and critical 
branes for A = fi. 

It is important to remark that the factor e~ A{ - r ^ (see eq. f|3T|) ) diverges for r — > ±oo in 
the case of DBW, such that it does impose restrictive conditions on the normalization of 
the wave functions. On the other hand, in the case of CBW, the factor e~ A ^ diverges 
for r — t- — oo and is constant for r — > oo (we have taken the upper signs in expressions 
(|T2l) -(( fl4l) ). Thus, the normalization of the wave function can only be found by means of a 
fine tuning on the Ui 's. 

In the case of CBW with A = fi one can see from (fT2|) and ( TIB"]) that the coupling 
W3 4>x contributes to a constant for the behavior of — 2r\ J r F{r')dr' for r — >• ±oo . In fact, 
we have found that 

t o r t?( >\a <\ J exp(-2r]acu 2 r), r +oo . . 

exp(-2?7 / Fir )dr ) ~ < ; _ ; , (45) 

J y exp(— zrjaujir), r — > — oo 

whilst e~ A ^ ~ e 4a3 Ml r l/ 9 f or r — oo. Then by choosing rj > 0, UJ2 > and u\ < —2a 2 [i/9r] 
one has localized left-handed massless fermion. On Figure 4 we show the profile for the 
r-dependent part of the coefficient of the spinor for lo\ = —1/3,-1,-3 and = 0, 77 = 
U2 = a = fj, = 1. As a conclusion, the coupling of fermions with the field 0(r) is relevant to 
provide localized massless fermions, but oj\ should be close to —2a 2 fi/9r] to insure a sharp 
localization in the core of the wall, since the peak of (r) dislocates to the right of the 
core of the wall as \u±\ increases. Moreover, the dominant contribution for r — > +00 comes 
from the coupling to the field x( r )> which describes the internal structure of the brane, 
and we have find that the coupling constant U2 must be positive in order to insure the 
normalization of the wave function associated to the massless fermion. We have also checked 
that the coupling u 3 4>x contributes to dislocate the peak of (r) from the core of the wall 
and that this coupling does not insure the fermion localization by itself, as has already been 
shown in [35]. Our results are in complete agreement with those found in reference [3T] . 



9 



In the case of CBW with A = 4/x one has the following asymptotic behaviors 



/r 
F(r')dr') 



exp(— 4rjau}2^), r — > +00 
exp(— 4r)auir), r — > —00 



_ AM \ const., r ->• +00 

e ( } ~ i / 16a • 46 

[ exp(- 9 M ), r -> -00 v ; 

such that one can choose 77 > 0, W2 > and Wi < —Aa 2 ^,/9f] in order to have normalizable 
fermion zero modes trapped inside the wall. 

We have also considered the localization of massless fermions on DBW with A = 4/x by 
taking the solutions (Q and fflUj) . We have found that smooth normalizable solutions are 
obtained by setting uj 2 = = and lo\ > 8a 2 ///97/. We show in Figure 5 the behavior of 
otL ( r ) f° r two different values of cq, Co <C 1/1 6a 2 and Co ~ l/16a 2 . One can note that the 
behavior of (r) does not follow the behavior of the brane. Whilst Figure 3 shows the 
two-kink solutions, which indicates the brane splitting, Figure 5 shows that the left handed 
massless fermion is likely to be seen in the bulk between the two walls as Co approaches 
the critical value. Although this seems to be a democracy in the sense that there is no 
preferable wall for the fermion to be trapped in, it also brings a paradox since the fermion 
lives in the bulk between the two walls preventing it to be measured. We have been looking 
for a solution to this apparent paradox. We have found that a coupling, reminiscent from 
supersymmetry, namely 

F(<t>, X) = ^1 TI/ — = wi — , (47) 

<p' 

with u>i < —2a 2 /9r] can afford localized massless fermion with definite chirality whose r- 
dependent behavior follows the brane splitting by exhibiting a sharp localization in the cores 
of the walls as can be seen in Figure 6. In this scenario the massless fermion can be found in 
both walls simultaneously and the probability density for the fermion to be found between 
both walls diminishes as the walls are far apart from each other. 



4 Conclusions 

In this work we analyze the localization of massless fermions on degenerate and critical Bloch 
branes. The calculations are done directly on the original extra-dimensional spatial variable, 
that is, the physical coordinates, in contrast to some works that resort to the transformation 
of variables dz = e~ A ^dr, that is, the conformal coordinates, in order to obtain effective 
time-independent Schrodinger equations, for the ^-dependent parts of the spinors, with eigen- 
value equals to the fermion mass. With such a procedure the effective Schrodinger equation 
allows us to see the emergence of a possible Kaluza-Klein (KK) tower of massive fermions. 
On the other hand, the appearance of possible non-localized KK modes is not evident from 
the effective Schrodinger equations in the physical coordinate r, equations (l4T)l) - (j4"3l . This 
is due to the fact that the mass contributes to the effective potentials, i.e., for each mass 
one has a different potential, such that one should look for scattering states, or even massive 
fermion localized states, by means of an effective Schrodinger equation with eigenvalue zero. 
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Notwithstanding, one does not have to resort to an effective Schrodinger equation in order 
to analyze the existence of normalizable localized states of massless fermions. 

It is important to remark that the interaction term of fermions with the scalar field is 
crucial to the correct localization of fermions on the brane. The most natural coupling is the 
Yukawa one, namely, (frtyty, used originally in to explain the charge fractionization by a 
soliton background, and stablished in [lj to illustrate localization of fermions in a domain- 
wall. In fact, the Yukawa coupling is the simplest one, but such a choice comes naturally 
if one has in mind the potential |(0 2 — l) 2 , whose BPS solution is the soliton background 
that traps the fermions. In this case, the Yukawa coupling entails a N = 1 supersymmetry 
(SUSY) in the fermion-boson system, once the superpotential is W(4>) = \/A0(^ — 1) and 
d 2 W/d(f) 2 = 2\/~\(j). In the context of branes in a warped space-time the issue of SUSY, that 
is supergravity, is much more complicate, moreover when it comes with two scalar fields, 
which is the case of Bloch branes in a model which supports a variety of soliton solutions. 
We have shown that the general coupling (ui(p + u 2 X + u 3 <fix)^^ guarantees the localization 
of massless fermions, for a range of the coupling constants ui\, u>2 and 0J3, and that the cross 
term (px does not provide fermion localized states by itself, in agreement with previously 
reported calculations [31], [35]. The general coupling also works even when one has two- kink 
solutions, but the r-dependent part of the wave-fuction is peaked just in the middle of the 
region between the walls. Such a behavior is not the desirable one, once it signalizes that 
the fermions would not be observed inside the branes. We have proposed another coupling 
between fermion and scalar fields which provides the correct localization of fermionic zero- 
modes inside the branes. The chosen coupling seems to come as a reminiscent of a SUSY 
model. This last issue has been analyzed in [36] . 
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,2 A (r) 




-1 



Figure 2: Warp factor in the case CBW with A = Afj, for /i — 1 and a = 1 (dashed line), 
a = 1.2 (thin solid line), a = 1.4 (thick solid line). 
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Figure 3: Two-kink solution (DBW) with A = 4//: c = (dashed line), c = 1/16.001 (thin 
solid line), c = 1/16.0001 (thick solid line) with a — ji — 1. 




Figure 4: «Lo( r ) in the case of CBW with A = 4/i and W2 = 1, 0J3 = 0, u)\ = — 1/3 (dashed 
line), ui\ — — 1 (thin solid line) and wi = —3 (thick solid line); a — /i — 1. 
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Figure 5: oloO") in the case of DBW (A = 4/i) and F = u\(p for ui\ = 4, c = (dashed 
line), c = 1/16.01 (thin solid line)and c = 1/16.0001 (thick solid line); a — ji — rj — 1. 




16 



